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The quantum statistics of damped higher-order optical solitons are ana- 
lyzed numerically, using cumulant-expansion techniques in Gaussian approxi- 
Q"^ ■ mation. A detailed analysis of nonclassical properties in both the time and the 

On . frequency domain is given, with special emphasis on the role of absorption. 

Highly nonclassical broadband spectral correlation is predicted. 



From classical optics it is well known that nonlinearities can compensate for the 
dispersion-assisted pulse spreading [^|| or for diffraction-assisted beam broadening (see, 
e.g., P). In the two cases, the undamped motion of the (slowly varying) bosonic field 
variables d(a;, t) is governed by the Hamiltonian 

g ; H = njdx {dxd^){dxd) + ^xa^a^dd] , (1) 



d{x,t),d^{x',t)] = S{x-x') (2) 



O 

GN ' [t, propagation variable; x, "transverse" coordinate; uj^'^\ second order dispersion or diffrac- 

tion constant; x nonlinearity constant; see, e.g., 0,0]. Note that bright temporal solitons 
can be formed either in focusing media with anomalous dispersion (x< 0, 0) or in 

defocusing media with normal dispersion {x > 0, 0), whereas spatial solitons require 

Ch ' always focusing nonlinearity. The effect of absorption is described in terms of ordinary 
^ ' Markovian relaxation theory resulting, in the low temperature limit, in the master equation 

ih dtp = [H, p] + i'jh J dx (2dpd^ — pa^d — d^dp^ (3) 

^ I (7, damping constant). 

■ The master equation (|^) is converted, after spatial discretization, into a pseudo-Fokker- 

Planck equation for an s-parametrized multi-dimensional phase-space function, which is 
solved numerically using cumulant expansion in Gaussian approximation [Q. The initial 
condition is realized by a multimode coherent state without internal entanglement, and it 
is assumed that the field expectation value corresponds to the classical A^-soliton solution, 
{d{x,to)) = N tto sech(a;/a;o), = 1,2,... (oq and Xq, mean amplitude and width of the 
fundamental soliton, respectively). 

Spectral properties can be studied introducing the Fourier-component operators 

1 

d(cj,t) = (27r)-2 / dxe"^''d{x,t). (4) 



Here we restrict our attention to correlations of photon number fluctuations. In the case 
of fiber soliton pulses the correlations in the cu-domain can be measured using appropriate 
spectral filtering (see, e.g., in |^). In the case of spatial solitonic beams the correlations in 
both the X and tu-domains, respectively, can be measured by filtering the field in the near- 
and far- field zones of the output beam (see, e.g., 0]). 
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FIG. 1. The evolution of the mean photon number (rij) and the correlation coefficient rja 
of an undamped soliton, N = 2, is plotted in the x-domain, Ax = 0.05 xq [(a), (6)], and the 
w-domain, ALi; = 0.25a;o [(c), (d)]. The plots (ei) — (es) (x-domain) and (/i) — (/s) (cj-domain) 
show the correlation coefficient ijij for typical propagation lengths {uq = I/xq, td = l^^o/'^^^^l' 
/ dxat(x,0)a(x,0) =8 x 10^). 



The output can be given by (see, e.g., |]^) 

biu, t) = G(z/, t)a(z/, t) + ^l-\Giu,t)\^ />, t), (5) 

where, according to the domain considered, u stands for x or u, and G{v,t), \G{v,t) \ < 1, is 
the (complex) transmittance of the fiher and /(z^, t) is a bosonic noise operator. The photon 
number operator of the detected hght is ri = / (iz/6^(z/, t)6(z/, t). Assuming square bandpass 
filters with Gi{v, t) = — Vti \ < Ai7 and Gi{v, t) = otherwise, we consider the correlation 

coefficient 




FIG. 2. The maximal violation of the Cauchy-Schwarz inequality for the photon number 
fluctuation [plots (a) — (d)] and the smallest Fano factor F= (An^)/(n) (strongest photon 
number squeezing) achievable with optimized filters [plots (e) — (h)] are shown for the funda- 
mental soliton, = 1, (dotted line) and the soliton with N = 2 (full line) [x-domain: plots 
(a), (6), (e), (/), w-domain: plots (c), (d), (g), (h); 7 = 0: plots (a), (c), (e), (5), 7id = 0.03: plots 
(6), (d), (/), {h); other parameters as in Fig. ||. 
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{rrii = (fii), Cij = (: AhiAfij :), Afii = ni — rrii), where : : introduces normal ordering. It can 
be shown that rju <1, and \r]ij\ < 1 for nonoverlapping intervals. A negative sign of the 
coefficient rju or a value smaller than unity of the Fano factor Fi = (Anf)/ {hi) = (1 — 
indicates photon number squeezing of the filtered light. 

From Fig. |I| it is seen that typical changes in the evolution of {hi) [Figs. |I|(a), (c)] and 
those of Tjii [Figs. |l](6), (rf)] and r]ij [Figs. ||(ei)— (cs), (/i)— (/s)] are closely related to each other. 
Near the points of soliton compression [maxima of ijii) in Fig. |I](a)] the formation of strong- 
correlation patterns is observed [Figs. ^(62) — (64), (/2), (/4)]. In contrast to the x-domain 
[Fig. |3(&)], sub-Poissonian statistics is observed in the tu-domain [Fig. |l]((i)]. Moreover, the 
correlation in the tu-domain extends over a larger interval (relative to the corresponding 
initial pulse width) than the correlation in the x-domain. One possible explanation of such 
strong, almost perfect correlation {\riij \ 1) can be seen in the instability of the classical 
A^-soliton solution. From a linearization approach the internal noise of a quantum 
soliton should be associated with interferences ||^ between the soliton components and the 
continuum part of the solution to the classical nonlinear Schrodinger equation, as obtained 
by means of inverse scattering method (see, e.g., |T^). The qualitative changes observed for 
turning from the fundamental soliton to higher-order solitons (A^ = 1 ^ A^ = 2, 3, . . .) are due 
to the presence of more than one soliton component. Discrepancies between the parameters 



(amplitude, group velocity, etc.) of the soliton components of the A^-soliton solution play 
the central role in establishing very strong internal correlations. 

Nonclassical correlation can be detected, e.g., by testing the Cauchy-Schwarz inequality 
for the normally ordered photon number variances. When it is violated, i.e., 

CiiCjj - 4 < 0, (7) 

then the photon number noise in the intervals i and j is nonclassically correlated. Figures 
^(a) — (d) reveal that the nonclassical correlation of the 2-soliton is substantially stronger 
than that of the fundamental soliton even for an absorbing fiber. Such an increase cannot 
be explained by a simple intensity scaling. The effect is obviously related to the mentioned 
instability of higher-order solitons. It is remarkable that there exist propagation distances 
for which the nonclassical correlation is stronger for an absorbing fiber than a nonabsorbing 
one. 

The strongest photon number squeezing (smallest Fano factor) achievable with an op- 
timized broadband filter is illustrated in Figs. 0(e) — (/i). Compared with the fundamental 
soliton, only a small increase of the effect is observed for the 2-soliton in the cu-domain 
[Fig. |^(e), 6.6 8.4dB]. On the contrary, a rather strong increase of the effect can be ob- 
served in the a;-domain [Fig. ^g), 3.3 — > 9.6dB], provided that losses can be disregarded. It 
is worth noting that the best photon number squeezing is achieved in the cu-domain for the 
fundamental soliton and in x-domain for the 2-soliton. The results show that the degree 
of squeezing sensitively depends on the domain considered. Hence, replacing the Fourier 
transformation in Eq. (^) [including Eq. (|])] with more general transformation that relates 
the fields in the two domains, may offer possibilities of further optimization. In particular, 
when we restrict our attention to linear transformations which can be realized experimen- 
tally by passive linear optical elements, then we are left with a two-dimensional integral 
kernel function to be optimized. In this way we may hope that also for other nonlinear 
quantum objects a considerable improvement of nonclassical features can be achieved. 
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